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Abstract 

The limit relations for the partial derivatives of the two-electron atomic wave functions at the two- 
particle coalescence lines have been obtained numerically using accurate CFHHM wave functions. 
The asymptotic solutions of the proper two-electron Schrodinger equation have been derived for 
both electron-nucleus and electron-electron coalescence. It is shown that the solutions for the 
electron-nucleus coalescence correspond to the ground and singly excited bound states, including 
triplet ones. The proper solutions at small distances R from the triple coalescence point were 
presented as the second order expansion on R and InR. The vanishing of the Fock's logarithmic 
terms at the electron-nucleus coalescence line was revealed in the frame of this expansion, unlike 
the case of electron-electron coalescence. On the basis of the obtained boundary solutions the 
approximate wave function corresponding to both coalescence lines have been proposed in the 
two-exponential form with no variational parameters. 

PACS numbers: 31.15.Ja, 31.15.-p, 31.10.+Z 
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I. INTRODUCTION 



Two-electron systems present an excellent basis both for testing the new quantum cal- 
culational approaches to many-body systems and for studying a number of photo-electron 
and other atomic processes. This is because such systems are the simplest ones with enough 
complexity to contain the main features of a many-body theory. This complexity arises 
from the electron-electron Coulomb potential which depends on the interelectronic distance 
^"12 =1 i"i — i'2 I- The proper Hamiltonian for infinite nuclear mass and charge Z, in atomic 
units used throughout the paper, is given by 



It does not depend on any experimental constants whose values change considerably with 
improvements in measurement equipment. Therefore, it provides a standard for theoretical 
calibration. 

It has been established Q| that relativistic and quantum-electrodynamic corrections to the 
energy levels of an atomic or molecular system require highly accurate nonrelativistic wave 
functions. Rayleigh-Ritz variational calculations provide a wave function with relative error 
approximately proportional to the square root of the relative error in the energy. Therefore, 
if the energies are used to estimate the quality of the wave functions, then it is necessary 
to calculate the nonrelativistic energies to far greater accuracy than would otherwise be 
needed. 

The alternative way for obtaining the very accurate and locally correct wave functions is 
a direct solution of the three-body Schrodinger equation. The Correlation Function Hyper- 
spherical Harmonic Method (CFHHM), employed in this paper, realizes successfully this way 
of solution P]. Accuracy of the method is comparable to the most sophisticated variational 
calculations. 

For problems in atomic or molecular physics, eigenfunctions of the Hamiltonian (1) exhibit 
singular behaviour at those points in configuration space where two or more charged particles 
come together and the resulting potential becomes infinite. For systems with zero total 
orbital momentum (S'-states) the wave function depends only on three scalar variables ri, 
r2 and ri2, i. e., \1/ = ^'(ri, r2, ri2). At the two-particle coalescences, the derivatives of the 
wave function \1/ have discontinuities characterized by the famous Kato cusp conditions j^, 




(1) 
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which have the simplest form for the S'-state of a two-electron atomic system : 

(9* 



dri 
dr2 



\r,=o=-Z^{0,R,R), [r2 = n2 = R] (2) 



\r,=o=-Z^iR,0,R), [n = n2 = R] (3) 



9^ 1 

^lm=o=2^(^'^'0). [n=r2 = R] (4) 

The conditions (2) and (3) pertain to the situation when the coordinates of one of the 
electrons and the nucleus coinside. These conditions are valid for the electrons, which have 
the same (triplet states) or the opposite (singlet states) spin directions. The condition (4) 
deals with coincidence of coordinates of two electrons. It is valid only for the singlet states, 
while due to the Pauli exclusion principle \l/(-R, R,0) = for the triplet states. The inclusion 
of functions with such cusps into trial wave functions has been shown to improve dramatically 
the rates of convergence of Rayleigh-Ritz variational calculations [J]. The using of the proper 
correlation function, which obey the Kato conditions (2)-(4), accelerates considerably the 
convergence of CFHHM approach, as well. 

It is known that the cusp conditions (2)- (4) are consequences of the Coulomb singularity in 
the potential and provide specific relations between the wave function and its first derivative 
at the points of coalescence. It was shown in Ref. [5] that the coalescence behaviour also 
uniquely determines the third derivative of the spherically averaged wave function in terms 
of the lower derivatives. The deduced relations are valid for any atom, molecule, or electron 
gas in any smooth external field. 

There are also singularities involving more than two particle, such as the triple-coincidence 
singularity in the Helium atom, when both electrons simultaneously approach the nucleus. A 
formal expansion in powers of the hyperradius r = ^/ rf + and its logarithm In r about this 
singular point was proposed by Fock 0] for the S'-state wave functions of the Helium atom. 
Subsequently, much efforts has been devoted to understanding this expansion. The 0(r°), 
O(r^) , and 0(r^ In r) terms in Pock's expansion are easy to obtain analytically. The O(r^) 
term in the expansion has been obtained in closed form by Maslen and co-workers, through 
the extensive use of computer algebra j?,!^. Myers and co-authors have examined their 
results, and have verified that the inclusion of this term in the expansion yields a continuous 



3 



"local" energy, whereas the "local" energy is finite but discontinuous at r = if the term 
is omitted. Forrey performed variational calculations of the ground-state energy of 
Helium. His basis set included an explicit treatment of the Fock expansion in hyperspherical 
coordinates and extra products of Laguerre polynomials in perimetric coordinates. This work 
has demonstrated that the use of Fock basis functions provided a substantial improvement 
in the convergence rate. 

We would like to emphasize that the calculation of the accurate wave function at the 
coalescence lines is a very difficult problem just because of their cusp nature. On the 
other hand, a number of atomic physics problems could be solved by using the functions 
appearing on the RHS of Eqs. (2)- (4). The processes of photoionization in the Helium atom 
and heliumlike ions |11] could serve as an example. 

It is well-known (see, e. g., 

BBS) 

, that using Hamiltonian (1) we can present the 
Schrodinger equation for ^-states of two-electron atom/ions in the form 



52^ ^52^ 2 5^ 2 9^ 4 9^ 
drf dr^ 9r^2 ^1 "^ri r2 dr2 dri2 



n r2 ri2 



riri2 J dridri2 \ r2ri2 J dr2dri2_ 

In this paper we provide the accurate analytic solutions of the Schrodinger equation (5) at 
the coalescence lines for both small and very large R. The Kato cusp conditions (2)- (4) are 
employed to solve the problem. 

II. ELECTRON-NUCLEUS COALESCENCE 

To investigate the case of coalescence of one electron and the nucleus in two-electron 
atom/ions, one should find the limit as, e.g., r2 approaches zero for both sides of Eq. (5). It 
is easier to perform this mathematical operation with the help of following relations: 

rl-r^ + r\2 „ r| - + r\2 



cosOi; -^—^ i^ = cos^2, (6) 

2riri2 2r2ri2 

where di is the angle between the vectors ri and ri2, and 62 is the angle between r2 and ru 
(see Fig. 1). It is clear that: 

lim 9i = 0; lim 62 = 7i/2. (7) 

1-2^0 '•2-»0 
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FIG. 1: Interparticle coordinates and associated angles for the two-electron atom/ions. 



Then, using Eqs. (6), (7), we can rewrite Eq. (5), taking the hmit as r2 approaches zero: 



+ 2 — ^ + 



drl 

— + 



drfo dridri2 



1 



+ 2- 



"■2=" R\dri dr 



12 



+ 



- + e]^{r,o,r). 



(8) 



Here we took into consideration that ri = ri2 = R as r2 approaches zero. Inserting the Kato 
condition (3) into the RHS of Eq. (8), and assuming that the terms of Eq. (8) must be finite 
in the whole 3-D space, we obtain 



52^ /^2^ ^2^ 

+ ^ + 2 7r^ + 



drf 



c?r?2 ' dridru 



2 (9^ 



1 - z 

R 



E ^(i?,0,i?). 



(9) 



We could certainly obtain the same equation if took the limit as ri approaches zero (r2 
ri2 = R). 
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It is seen that the LHS of Eq. (9) presents a sum of the form J2kim ^kim(R)^im(-^)i where 
for the partial derivatives of ^ taken at the electron-nucleus coalescence line, we have 

^^^^^ = !^o^-^^^^ ^'^ = 0'l'2.12). (10) 

Here dro — 1. Then, in the case of the first partial derivatives we have k — 1, whereas one 
of the numbers Z, m is equal to zero. 

Let us now denote the two-electron wave yS-function at the electron-nucleus coalescence 
hne as 

^{R,0,R) = F{R). (11) 

So, if we could express all of the functions Af'liR) through the functions F{R), F'{R) and 
F"{R), with factors being depending on R, then we obtain an ordinary differential equation 
of the second order for the function F{R). The prime denotes differentiation, as usual. 

Solution of this differential equation under the proper boundary conditions could give us 
the desired function F{R). We do not yet know how to do this in general form. However, 
as a first but important step we propose here a method for solving Eq. (9) in the boundary 
regions, i.e., at very large R and at small R. One should notice that the numerical calculation 
of F(i?) in these regions is particularly difficult. 

The direct Correlation Function Hypcrspherical Harmonic Method allows to calculate 
numerically the two-electron wave function \E'(ri^r2,ri2), as well as its special case F{R), 
with very large accuracy. By using the CFHHM numerical calculations we obtained the 
following limit relations between the functions mentioned above for the asymptotic region 
of very large R: 

= 1, (12) 
= 1, (13) 



(14) 




0, (15) 
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Ar2!i2(^) 

F{R) 



FiR) 




(16) 



(17) 



The calculations show that these relationships arc valid at least to four significant digits. We 
cannot achieve higher accuracy due to the fact that the inaccuracies of the wave functions 
and especially of their derivatives go up with R. Note that the asymptotic relations (12)- (17) 
are valid for the ground states of the two-electron atom/ions, as well as for its excited states, 
including triplet states. Of course, we were not able to perform the proper calculations for 
all of the excited states. However, we could guarantee the validity of Eqs. (12)-(17) for 
several first ones. 

Relations (15)-(17) show that we can neglect the partial derivatives with respect to r^. 
Moreover, the calculations of the accurate CFHHM wave functions show that the ratio of 
F'{R)/F{R) achieves a finite value as R approaches infinity. In the next sections we will 
obtain this ratio as a finite function of Z and E. This property together with the limit 
relation (12) allows us to neglect the terms proportional to 1/i? on both sides of Eq. (9). 
And finally, using the relations (13), (14), Eq. (9) is transformed into the following simple 
differential equation: 



The function Fas{R) is the asymptotic representation (for very large R) of the accurate two- 
electron wave function in the situation when one of the electrons "seats" on the nucleus, 
while the other electron is far away. Eq. (12) shows that this function depends on two 
parameters, the nuclear charge Z and the total energy E of the two-electron atomic system. 
Ci is an arbitrary constant. We used only the accurate wave functions of the discrete 
spectrum {E < 0) to obtain the relations (12)-(17). The condition of exponent in Eq. (19) 
to be real leads to the inequality 



at least for the ^'-states of the Helium atom or hehumlike ions with the nuclear charge Z. 




(18) 




(19) 



-2E > Z^ 



(20) 
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We have calculated the CFHHM energy levels of the Hehum atom and the ions of Li+and 
for both the singlet and triplet ^'-states with n <7. These data, presented in Table I, 
confirm the validity of the inequality (20). 

TABLE I: Energy levels (in a. u.) of n^S states for the Helium atom and heliumlike ions of 

Li+{Z = 3) and B^+{Z = 5). 



n \ 


He: k = l k = 3 


Li+ : k = l k = 3 


: k = l k = 3 


1 


-2.903721 


-7.27887G 


-22.02788 


2 


-2.145970 -2.175225 


-5.040179 -5.110019 


-14.57652 -14.73188 


3 


-2.061221 -2.068696 


-4.733102 -4.751430 


-13.41017 -13.45127 


4 


-2.033566 -2.036524 


-4.629208 -4.636571 


-13.00797 -13.02461 


5 


-2.021225 -2.022633 


-4.581895 -4.585572 


-12.82326 -12.83163 


6 


-2.014537 -2.015122 


-4.556331 -4.558559 


-12.72341 -12.72824 


7 


-2.010629 -2.010870 


-4.541111 -4.542445 


-12.66341 -12.66654 



It is easy to conclude that Eq. (20) corresponds to the first ionization threshold Ii — Z^/2, 
and that the limit relations (12)-(17) describe the ground and singly excited bound states 
[12,13]. So, if the electrons are far away from each other (R oo), then the simplest model 
one may think of is the model of two independent electrons. The inner electron is bound in 
a state with principal quantum number N and energy En — —Z'^/{2N'^), the outer electron 
is in a hydrogen-like orbital with energy — —{Z — l)^/(2n^) and n > N assuming a 
screening of the nuclear charge. The total energy is simply the sum of the one-particle 
energies. The case of R approaching infinity corresponds to the conditions of n oo, 
En 0, and Ethresh = h = En=i- So, these arguments give additional evidence to the 
validity of the limit relations (12)-(17) and, consequently, the asymptotic solution (19) for 
the ground and singly excited bound states. 

We have obtained some limit relations for the functions (10) and (11) in the vicinity of 
the triple collision point, R — 0, ets well. The proper numerical calculations yield for the 
singlet states: 

, JaS>)1 ^ 
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Ita^^Ui, (22) 

^ F{R) J 2' ^ ^ 

Note that Eq. (21) contains the singlet function F{R), whereas the corresponding Eq. (12) 
includes the first derivative F'{R). The results for the triplet states are: 

= (26) 




F(i?) I 2' 
Ago(^) 1 1 



(27) 



F{R) J 4' 

As is known, the two-electron wave functions of the singlet states at the triple coincidence 
point are non-zero. Therefore, Eqs. (21), (22) allows to avoid the divergence at i? = for 
terms proportional to R~^ in Eq. (9). The triplet states, which are proportional to R^ in 
the vicinity of the triple collision point, don't have to obey such a requirement. 

We were able to obtain only the simplest limit relations as R approaches zero. We hope 
these relations could be good for searching the general solution of Eq. (9). However, to 
derive the solution for small R we propose another way, which is more precise as well as 



analytic expansions of the three-body 



more reliable. As was mentioned earlier, in Refs. 
atomic wavefunctions were presented. The expansions were derived for the exact solutions 
of the Schrodinger equation (5) up to the terms of the order (including Inr). Note that 
limit relations (21)-(27) {R 0) could be obtained by using Ref. [8] too. We used some 
of those results (see Ref. [8], pp. 2796-2797) to obtain the analytical representation of the 
two-electron wave functions at the two-particle coalescence lines in the vicinity of the triple 
coincidence point. The same results, but for the singlet states only, could be derived by 
using Ref. [9]. However, one should be very careful, because in the last reference we found 
at least three misprints, which could have influence on the final results. The first misprint 



is a missing factor 2 in the expression for 1^2,0 (below Eq. (14) j^). The second one is the 
incorrect expression ri2 sin a cos on the RHS of the expression for 1^2,1 (below Eq. (14) j^). 
And the third misprint is the missing function cos~^ before (ri ■ r2/rir2) in the RHS of Eq. 
(11) 3. 

So, using the results of Refs. [8,9] and taking the limit as r2 0, we obtain for the singlet 
states: 

/ 1\ i?2 

F{R) ~ 1 - i? f Z - - j + Y^[4^^ - 2Z(3-ln2) + 1 ~ 2E]. (28) 
The similar result for the triplet states has the form 



2„ 1\ H 



?2 



-Z2-Z(2--ln2) + --E 
3 ^6^4 



(29) 



For simplicity, the wave functions (28), (29) are normalized by condition -F(O) = 1 for the 
singlet states, and [F(_R)//2^]^^q = 1 for the triplet states. 

We would like to pay particular attention to the absence of the Fock's logarithmic term 
in both expressions. This term disappears at the electron-nucleus coalescence line, because 
of the vanishing Fock's angular coefficient il'2i in the limit as r2 approaches zero (whereas 
ri ^12) for the singlet states. 

At first glance, it is natural to assume that all of the Fock's logarithmic terms are canceled 
at the electron- nucleus coalescence line. However, such an assumption proved to be incorrect. 
We verified the angular coefficients ip^i, ipn and ■?/'42, corresponding to the logarithmic terms 
in the Fock's expansion up to the terms of order r'^, r'^ln^r (singlet states). The exact 
expressions for these quantities could be found, e.g., in Ref. [7]. All of these three angular 
coefficients proved to be nonzero. So, we conclude that the first logarithmic term il)2ir'^\n.r 
of the Fock's expansion is the only one to vanish at the electron-nucleus coalescence line, 
at least for the singlet states. Eq. (29) shows additionally that for the triplet states all 
of the logarithmic terms, up to the fourth order in r, are canceled in the limit as r2 ^ 
(or ri — > 0). Accordingly, the values of the first and second derivatives F'{{]), F"(0) for 
the singlet states and [F(_R)/_R^]'^^q, [F(_R)/_R^]'^^q for the triplet states are finite. We 
verified the validity of the expansions (28), (29) by direct calculation of these derivatives 
in the limit as i? — > 0, using the accurate \l/-functions. The results coincided with the 
calculations performed according to the analytical formulas (28), (29) within the accuracy 
of five significant digits. 
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III. ELECTRON-ELECTRON COALESCENCE 



For the case of forming the two-electron coalescence or the coincidence of the coordinates 
of the both electrons, one should take the hmit as ri2 — > on both sides of the Eq. (5). In 
this case we have (see Fig. 1): 

lim 9i = lim 62^7^/2. (30) 

'"12— >0 T-12— >0 

Then, in the limit as ri2 — > both terms with the mixed partial derivatives vanish in Eq. 
(5), and we can write: 

1 / , 2 d-^ 

2^ + 2- 



2\ drl drlJ ''■^'=° Rdn 

-M L.^. + (^4 + E)m,R,0). (31) 



ri2 \dri2 2 J ''-1^=° V R 

When deriving Eq. (31), we took into account that ri = r2 = i? as approaches zero. 
Then, using the Kato condition (4) in the RHS of Eq. (31), we obtain 

+ ^ In.^o +-5^ U.o= - -^ + e)^{R,R,0). (32) 



drf drj^J Rdn ''^^ \R 

From the CFHHM numerical calculations, we obtain the following relation 

y a^(ri,r2,ri2) 1 d^R) 

^™ a = o ,p > 33) 

or I 2 dR 

where 

$(i?) = *(i?,i?,0). (34) 

The relation (33) is not a double-limit relation like all the other similar relations presented 
in the Section II. It is the only relation of such form obtained. It is valid for all < i? < 00 
and Z >1. The excited states are included as well. Using Eq. (33), we can rewrite Eq. (32) 
in the form 

52^ 
drf dr^2 



2 +^]U.^.= 9{R). (35) 



where 

Calculation of the ordinary first derivative d^/dR is considerably more precise than the 
calculation of the partial derivatives of the second order. The proper calculation of the 



l^^ + E^^R). (36) 
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function g{R) and its comparison with the ordinary second derivative d?(^/dR^ yields, in 
the hmit of R approaching infinity: 

lta|^;^l=4. (37) 



S(R) 

Then, using the hmit relation (37), we can rewrite Eq. (36) in the limit of very large R: 



- as 



+ E<^as = 0. (38) 



4 dR^ 

Like before, we neglected the terms proportional to R^^ in Eq. (36). The proper solution of 
the differential equation (38) has the form 

^as (R) = Ci exp (-2/?v/=e) , (39) 

where Ci is an arbitrary constant. The function ^as{R) is the asymptotic representation of 
the accurate two-electron atomic wave function for the case of two electrons being at the 
same point but far away from the nucleus. Unlike F^s, Eq. (39) shows that depends on 
only one parameter, energy E. As far as we consider the discrete spectrum, then E < and 
consequently the exponent on the RHS of Eq. (39) is negative. 

Moreover, we obtained the following relation in the limit as i? — > 0: 

lim lim = -Zlim$(i?). (40) 

It enables one to avoid the divergence at i? = for the terms that are proportional to R~^ 
in the general Eq. (32) for the case, when two electrons are at the same point. 

The approximate solution of the Eq. (32) at small R may now be obtained by neglecting 
the terms that are not proportional to R~^ in that equation. Using also Eq. (33), we have 
the equation 

S + 2^'^2 = 0, (41) 

with the the solution 

<^2{R) = C2exp{-2ZR). (42) 

We can again use the approach developed in the works 0, 0] at ri2 — (ri = r2 = R). 
However, one shouldn't forget that the expressions presented in these papers are valid only 
for r2 < ri. In Ref. the angular coefficients '^/'ocbj^io, "^21, "^20 of the Fock expansion are 
expressed through Pluvinage coordinates C and r] jl^, which are more symmetric and are 
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valid for arbitrary ri and T2- Taking into account that the scahng transformation r — r/Z 
was apphed to the Hamiltonian in Ref. [14], we obtained under the simphfying normahzing 
condition of $(0) = 1: 

~2Z (2 \ , „ E 5 



~ 1 - 2ZR + i? 
with 



1 1 Ini?- - + -Z^ + 4 
3 Vvr / 3 3' 



(43) 



1Z / I \ \ Z 

Here 02,1 is an unknown coefficient of the homogeneous solution, which cannot be deduced 



from the local behaviour of the Schrodinger equation near the nucleus ja |9|, |10[ . It is a very 
difficult problem to calculate exactly 02,15 because of the presence of the logarithmic term 
on the RHS of Eq. (43). However, we have evaluated 5z (and consequently 02,1) using the 
accurate CFHHM wave functions. This yielded the following approximate dependence: 

5, ~ - - —Z. (45) 
2 10 ^ ' 

This linear part of the Z-dependence is the most accurate one for 2 < Z < 5. We found 
that the term —Ej'i in (43) expresses correctly the dependence of $(-R) on the state of 
excitation, while the term (5/3)Z^ represents accurately the nonlinear dependence of $(-R) 
upon the nuclear charge Z . Note that using the proper formulas from Refs. [8,9] yields the 
erroneous term Z^ instead of the correct result (5/3) in the expression (43). As we can 
see, the Fock's logarithmic term is preserved at the electron-electron coalescence line, unlike 
the case of the electron- nucleus line. Hence, the second derivative $"(i?) has a logarithmic 
singularity at the origin. It is seen that solutions (42) and (43) coincide in the first order 
approximation at small R. 



IV. RESULTS AND DISCUSSIONS 



In Sections II, HI we have obtained analytical representations of the accurate two-electron 
wave functions at the boundary regions of both two-particle coalescence lines. We found 
that the behaviour of these boundary solutions both at small and at large distances R has 
an exponential character. These properties enable us to propose a simple approximate rep- 
resentation for the ground state wave functions in the two-exponential form. The main idea 
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is that the first exponential represents the behaviour at small R and the second exponential 
represents the behaviour at very large R. So, let us consider the following function: 



f{R) = C{exp(-Ai?) + 7exp(-/3i?)} , (46) 

with 

A = «(! + 7) - 7/3. (47) 

Here C, a, /5 and 7 are arbitrary parameters. This two-exponential function possesses an 
important peculiarity: 

that is the ratio of the first derivative to the function itself at the origin (R — 0) depends 

upon only one parameter a. Let the second exponential in the RHS of Eq. (46) present 
the behaviour of the wave functions in question at very large R. Then according to the 



results of the previous sections (see Eqs. (19), (28)), we can put Ps = V—'^E — Z"^ and 
Us — Z — 1/2 ior the electron- nucleus coalescence fine of the singlet states. Considering C 
as the normalization constant, we have only one unknown parameter, 7. It can be obtained 
using the second derivative of the wave function F{R) at the origin. Then the double 
differentiation of the general function (46) yields, in the limit of R approaching zero, 

m~ 1+7 ■ 

On the other hand, we have from Eq. (28): 

hs^^ [4^^ - 2Z(3 - In 2) + 1 - 2£;] . (50) 

Eq. (49) is a quadratic with respect to parameter 7. It has two roots: 71 = —1 and 
72 = (/i — a'^)/{a — py. While the root 71 yields a trivial solution f{R) — 0, we retain only 
the solution 7 = 72. Using Eq. (50) and taking into consideration the values of parameters 
as and fig mentioned above, we obtain 

A(Z \n2 - E - Z^) - 1 , , 

Is = -f , ' ^,2 (51) 

3 (1 - 2Z + 2V-2^ - Z2) 

The exponent A of the first exponential plays an important part in constructing the approx- 
imate function (46). Note that it does not represent the behaviour of this function at the 
origin. Both exponentials give contributions to the correct behaviour of the function (46) 
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in this boundary region according to Eqs. (48), (49). However, the inequahty \ > (5 must 
be vahd if we want the exponent (5 to present the approximate function (46) at very large 
R. So, using Eq. (51) and the proper expressions for parameters and (5s we have for the 
ground state 

_ 2 - 4£; + - 4Z[3(1 + (5s) - In 2] + 6/3, 

In Table II we present the numerical values of the exponents A, and (5s as well as the factor 
7s for the Helium atom and several heliumlike ions. As one can see from Table II, the 



TABLE II: The parameters A, (5 and 7 for the approximate wave functions at the electron-nucleus 
coalescence line. The lower index s signifies that the parameter presents a singlet state. 



z 


2 


3 


4 


5 


As 


1.58574 


2.75657 


3.93411 


5.11350 




1.34441 


2.35793 


3.36320 


4.36599 


7^ 


0.551062 


1.80594 


3.17337 


4.57790 



parameter A^ is greater than (5s for all Z presented. 

To verify how good are the approximate functions (46) with the parameters of Table 
II, we have drawn the curves y{R) = log^o ^"'^^^r)^^ Fig. 2. Fs{R) and fs{R) are 
the accurate and the approximate functions, respectively. One can see from Fig. 2 that the 
approximate curves are very close to the exact ones at small R. However, even at quite large 
R the accuracy is not lower than 10%. The total accuracy is increasing with the nuclear 
charge Z. 

The triplet states of the two-electron atomic systems are always excited states. Only the 
electron-nucleus coalescence line is formed for these states according to the Pauli exclusion 
principle. It is important that the corresponding wave function at the coalescence line -ft(-R) 
behaves according to Eq. (29) like R^ as R approaches zero. Therefore, the function (46) 
is not suitable in this case. Instead, for the lowest energy triplet state we can propose the 
simple approximate function of the form 

MR) = C,R[eM-\R) - exp(-Ai?)], (53) 
where (5t — \J —2Et — Z"^ has to describe the behaviour of F{R) at very large R. It is easy 
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FIG. 2: Deviation of the approximate function fs{R) from the exact value Fs{R) at the electron- 
nucleus coalescence line for the ground states of the two-electron atom/ions with the nuclear charge 
Z = 2,3,4,5. 



to derive the following property of this function at the origin 




ftjR) 

i?2 



X 



jm 

i?2 



(54) 



On the other hand, according to Eq. (29) the exact representation of Ft{R) yields the value 
of (1/4 — 2Z/3) for this quantity. So we obtain for the unknown exponent 



(55) 



It is easy to verify that the inequality At > is valid for all Z. This condition is necessary to 
make the behaviour of ft{R) close to the correct one. Fig. 3 demonstrates quite satisfactory 
behaviour of the approximate functions (53), especially for i? < 5 a.u. We can observe that 
at small R the accuracy of ft{R) decreases with increasing Z, whereas at large R it increases 
with Z. The dips on the graphs of Fig. 2 and 3 are artifacts of the logarithmic scale, since 
the logarithm of the absolute value of the difference of the two functions goes to — oo at the 
points of crossing the functions. The overall accuracy therefore can be inferred only at the 
values of R not too close to the dip. 
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FIG. 3: Deviation of the approximate function ft{R) from the exact value Ft[R) for the 2S'-triplet 
states of the two-electron atom/ions with the nucleus charge Z = 2,2,, 4, 5. 

The case of the electron-electron coalescence is the most complicated, because of the 
presence of the logarithmic term in the Eq. (43). However, the general two-exponential form 
(46) can be applied in this case too, but without the representation (47). 

For simphcity, let us put C = 1/(1-1-7) Eq. (46). Then we obtain the approximate 
function fe{R), which obeys the condition /e(0) = 1, that in turn enables one to compare 
fe{R) with the expansion (43) in a simple manner. The exponent j3 representing the be- 
haviour at very large R is equal to I3e = 1\J —E and has to obey the inequality Ae > /9e, 
according to the asymptotic representation (39). So, we have two unknown parameters Ae 
and 7e. The factor 7e has to be positive, because the ground state wave function is nodeless. 
Then, for the first derivative fJyRi) at the origin we have /e(0) = — (Ag + /3eAe)/(l + 7e)- 
Comparing this formula with Eq. (43), we obtain the first constraint: 



It is easy to verify that the inequality 2Z > (5^ is vaUd for ell Z > 1. Therefore, for the 
parameters Ae and 7e obeying the constraint (56), the condition Ae > (3e will be valid for 
any positive 7e. To obtain the second constraint for the parameters Ae and 7e, one can use 



Ae = 2Z + 7e(2Z-/3e). 



(56) 
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for example an integral property of the exact wave function $(-R) such as the normalization 
integral S = ^^{R)R^dR with $(0) = 1. It is easy to calculate S using the accurate 
CFHHM wave functions available. Replacing $(-R) by fe{R) in the integrand, and executing 
a simple integration, we obtain the second constraint in the form 



27e 



+ 



Te 



(57) 



The simplest way of solving the set of the equations (56), (57) is using the well-known 
program Mathematica Q]. The equations have a number of roots (including complex ones). 
However, only one root turned out to be real and positive, and therefore it could be applicable 
to 7e. The corresponding solutions along with the other accompanying quantities, E, S 
and are presented in Table III. It is seen from Table III that the approximate function 
parameters fulfill all the conditions mentioned above. 



TABLE III: Parameters /?e, Ae, and 7e of the approximate wave functions at the electron-electron 
coalescence line. The accompanying values of the total energy E and the integrals S are presented 
as well. 



z 


-E 


S 




Ae 


7e 


2 


2.9037244 


0.00452855 


3.40806 


5.54012 


2.60184 


3 


7.2799134 


0.00127463 


5.39626 


8.32976 


3.8589 


4 


13.655566 


0.000524535 


7.39069 


11.1232 


5.12583 


5 


22.030917 


0.000264651 


9.38742 


13.9178 


6.39552 



The degree of coincidence for the accurate and approximate Helium wave functions at 
the electron-electron coalescence line is presented in Fig. 4. The solid line describes the 
deviation curve for the two-exponential approximate function fe{R) of the Helium atom. For 
comparison, we have also drawn the corresponding one-exponential function (42), presented 
by the dotted line. The graphs are limited by the value of i? = 2 a.u., because of the very 
fast decay of the two-electron atomic wave functions at the electron-electron coalescence line. 
For example, if $(0) = 1, then $(2) ~ 0.0006 (for Z = 2). For comparison, we can point 
out that the corresponding value of the wave function at the electron-nucleus coalescence 
line Fs{2) ^ 100 $(2). The graphs on Fig. 4 demonstrate that the two-exponential function 
is considerably more accurate than the one-exponential one, $2(-R)- AH of the approximate 
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FIG. 4: Deviation of the approximate functions fe{R) and ^2{R) from the exact wave function 
at the electron-electron coalescence line for the Helium atom. 

functions presented in this Section could be employed for the evaluation of different atomic 
phenomena, and processes of the atomic photoionization in particular. 

V. CONCLUSION 

We have considered and analyzed some particular solutions of the Schrodinger equation 
for the two-electron atom or ion (with the nucleus charge Z and the total energy E) using 
the accurate CFHHM wave functions. We have obtained mathematical relations between the 
partial derivatives taken at the two-particle coalescence lines and the ordinary derivatives 
of the wave function taken at the same coalescence lines. The relations were found for the 
hmit cases of very large and small distances R between one of the electrons and the other 
electron close to the nucleus (electron-nucleus coalescence) or between the two electrons 
close together and the nucleus (electron-electron coalescence). We have obtained the only 
relation valid for all i? > and Z > 1. It connects the first partial derivative on ri (or r2) 
and the first ordinary derivative of wave function at the electron-electron coalescence line. 
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We have examined both singlet and triplet S'-states, and obtained the asymptotic so- 
lutions (for very large R) with the exponents [—R\/—2E — Z^^ or {—2R\/—E) for the 
electron-nucleus or electron-electron coalescence, respectively. These results turned out to 
be valid for the ground and singly excited bound states. Wc have derived the second order 
expansions in R and In R in the vicinity of the triple coalescence point for small R. We have 
found that the Fock's logarithmic terms vanished at the electron-nucleus coalescence line, 
unlike the case of the electron-electron coalescence, but only in the framework of the second 
order expansion. 

We have proposed simple two-exponential approximations for the two-electron atomic 

wave functions at the coalescence points. The approximations arc valid for the lowest energy 
(nodclcss) states. It was demonstrated that all of the approximate functions have quite 
satisfactory accuracy and could be employed at least for the initial qualitative evaluation of 
a number of phenomena in the atomic physics. 
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